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' Key words, galaxies: clusters: individual (Coma) 
1. Introduction 

Long considered the archetype of rich and relaxed clusters, 
Coma is now known to have a rather complex structure, 
with a number of substructures detected at various wave- 
, lengths (see Biviano 1998 for a complete review on the 
' Coma cluster up to 1995). Indeed, although a large quan- 
. tity of multi-wavelength observational data are presently 
' available for Coma, our theoretical understanding of this 
unique cluster is still far from complete. 

Previous photometric studies of the Coma cluster were 
done in few colors, and were limited either to wide but 
shallow surveys or to small deep ones, often obtained un- 
der less-than-ideal seeing conditions. For example, very 
large shallow surveys (e.g. Godwin & Peach 1977, here- 
after GP77) were the basis for important studies of the 
Coma cluster luminosity function and of the properties of 
galaxies in general. However, due to the shallowness of this 
catalog, very little was known about Coma's galaxy pop- 
ulation at magnitudes fainter than R~20. Recent deeper 
catalogs are still limited to R'-~^21 when covering the whole 

* Based on observations obtained at the Canada-France- 
Hawaii Telescope (CFHT) which is operated by the National 
Research Council of Canada, the Institut des Sciences de 
rUnivers of the Centre National de la Recherche Scientifique 
and the University of Hawaii. 



cluster (e.g. Castander et al. 2001, Terlevich et al. 2001 or 
Komiyama et al. 2002). Finally, very deep catalogs such 
as that of Bernstein et al.(1995) covered only a very small 
area of sky. While many of these previous studies led to 
important discoveries, all suffered to some extent from in- 
sufficient depth or coverage, which made it difficult to gen- 
eralize the results obtained from these studies. 

Here we report a new catalog of Coma cluster galax- 
ies, obtained with the CFH12K camera (Cuillandre et al. 
2000) installed at the prime focus of the Canada-France- 
Hawaii Telescope (CFHT), which provides the first deep, 
wide-field images of Coma in the B, V, R and I filters, in 
seeing conditions well below 1 arcsec for a large part of 
the observations. As the field covered by the CFH12K is 
0.32 deg^ we made two overlapping pointings to cover a 
total field of 0.6 deg^. The data obtained are significantly 
deeper than any of the recent CCD large field imaging of 
the Coma cluster; e.g. Terlevich et al. 2001 (U and V fil- 
ters, depth V=20, 1 deg^, 2.2 arcsec seeing) or Komiyama 
et al. 2002 (B and R filters, 2.25 deg^ complete to R=21, 
seeing from 0.8 to 1.5 arcsec). A BVR trichromic image of 
the final result is shown in Fig. ^ 

Our main scientific goals in obtaining deep wide-field 
multi-color images of Coma are: a) to study environmen- 
tal effects on the faint end slope of the galaxy luminosity 
function (as a very significant extension of, e.g., the work 
by Lobo et al. 1997 or by Andreon & Cuillandre 2002); 
b) to analyze Coma's morphological and dynamical struc- 
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Fig. 1. BVR color image of the Coma cluster from our CFH12K data. The size of the field is 0.72degx0.82deg. North 
is up and East is to the left. 



ture to see how it relates to the cluster's formation history 
(e.g. Melher et al. 1988 or Adami et al. 2005b), the dif- 
fuse light emission (Adami et al. 2005a) and the faint and 
low surface brightness galaxy distribution (Adami et al., 
2006, submitted). However, because this new Coma cata- 
log provides a rich source of data that will undoubtedly be 
useful for other researchers, we intend to make it publicly 
available. 

In this paper, we describe our data acquisition, 
reduction, and cataloguing. Catalogs are available at 



Table 1. Main characterictics of the detectors 



http://cencosw.oamp.fr/ 



2. Observations and data reduction 

2.1. Requirements 

Several considerations guided our survey design; we sum- 
marize these below. 

Requirement 1: One of our primary science drivers was 
to investigate environmental effects on the galaxy popula- 
tion in the Coma cluster. Limiting ourselves to the densest 
part of the cluster, we wanted to cover at least 5 core radii, 
in order to sample regions where the density falls below 
~1% of the central cluster galaxy density. To do this re- 
quires surveying a region 50x50 arcmin^, which is possible 
with two adjacent CFH12K fields (see next section). 

Requirement 2: We study the spectral energy distri- 
butions of galaxies in order to constrain galaxy formation 
models. Four photometric bands were required for this, in 
order to obtain color-color plots. We also wanted to com- 
pute photometric redshifts and this requires, at the cluster 
redshift, to have at least one band below the 4000 A break 
and one above it. For Coma, this means that we needed 
U band data (not yet acquired) and B, V, R and I band 
data (all observed). 

Requirement 3: We wanted to sample the dwarf galaxy 
population as deep as possible in at least one band, with- 
out suffering too much contamination from globular clus- 
ters. At the redshift of Coma, this implies needing a com- 
plete catalog down to R=24 (see Bernstein et al. 1995). 

2.2. Instrumental setup 

Our data were acquired at the CFHT 3.6m telescope 
with the CFH12K CCD mosaic (12 individual 2048x4096 
CCDs) (see Cuillandre et al. 2000 and Table [ij installed 
at the prime focus. The pixel size is 0.206 arcsec, well 
suited to the mean seeing at prime focus (^^0.8 arcsec), 
giving a full field of view of 42x28 arcmin^ per field. Two 
adjacent fields offset in the North/South direction with a 
'^7 arcmin overlap were observed in the B, V, R and I 
CFH12K filters giving a final field of view of 0.72x0.82 
deg^ (in good agreement with Requirement 1). The V, R 
and I are Mould filters while B is a CFH12K customized 



Detector 


CFH12K 


Pixel size 


0.206 arcsec 


Gain 


1.4-1.8 e-/ADU (typically 1.5) 


Quantum efficiency 


80% at [600,700]nm 


Readout noise 


3 to 10 e~ (typically 5) 


Dark current 


1 e~ per min 



Table 2. Filter characterictics 



Filter 


A(A) 


AA(A) 


Peak transmission 


B 


4310 


950 


85% 


V 


5370 


940 


90% 


R 


6580 


1300 


85% 


I 


8223 


2164 


91% 



filter (see Table We will refer hereafter to the filters 
using the capital letters B, V, R and I and to the fields 
as North for the North pointing and South for the South 
pointing. These four bands partially satisfy Requirement 
2, and we still plan to acquire the missing U band data. 

2.3. Exposure times and first/second epoch data 

A first set of B, V, R observations was obtained in 1999 
and 2000 with only 10 of the 12 CCDs available at that 
time (first epoch data). A second set of CFH12K images 
(second epoch data) with all CCDs working, was obtained 
in 2000: the South field only in R, both North and South 
fields in L In both runs, long exposures (1 to 2h) were split 
into several shorter ones, offset by a few arcsec, in order 
to improve the cosmetics, flat field and fringing correc- 
tions. Offsets between exposures were kept small enough 
so that images could be stacked without correcting for the 
differential distortion between exposures. 

Exposure times were 24 min in B, 12 to 14 min in V, 
10 to 12 min in R and 12 min in I, rather long (which 
resulted in bright galaxies being saturated) in order to be 
sky background limited and to minimize the overheads. 
Tables Bl Eland ^ summarize the main observational 
informations. 

2.4. Data reduction 

Acquiring and reducing the data was part of C. 
Savine's PhD Thesis (Savine 2002) and the process 
is fully described in her thesis report (available at 
http://cencosw.oamp. fr/l. Here we summarize the key 
steps. 



First epoch data were preprocessed by J.C. Cuillandre 
using his own FITS Large Images Processing Software 
(FLIPS) package, which performs functions similar to 
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Table 3. First epoch observations 



Fiher 


Field 


Center (J2000) 


Exposures 


Total exposure time 
(s) 


Mean seeing 
(arcsec) 


Observing vear 


B 


South 


12:59:40.00 


27:48:47.1 


5 


7200 


1.07 


2000 


B 


North 


12:59:40.00 


28:10:07.0 


5 


7200 


1.01 


2000 


V 


South 


12:59:40.00 


27:48:47.1 


5 


4200 


1.00 


2000 


V 


North 


12:59:40.00 


28:10:07.0 


7 


5040 


0.97 


2000 


R 


South 


12:59:40.00 


27:49:07.2 


6 


3300 


0.90 


1999 


R 


North 


12:59:40.00 


28:10:06.9 


5 


3600 


0.87 


2000 



Table 4. Second epoch observations 



Filter 


Field 


Center (J2000) 


Exposures 


Total exposure time 


Mean seeing 


Observing year 










(s) 


(arcsec) 




R 


South 


12:59:40.00 ; 27:49:17.0 


10 


7200 


0.85 


2000 


I 


South 


12:59:40.00 ; 27:48:46.9 


10 


7200 


0.80 


2000 


I 


North 


12:59:40.00 ; 28:10:06.9 


10 


7200 


1.00 


2000 



IRAF task MSCRED and is optimised to handle large 
mosaics. 

Second epoch data were fully preprocessed at 
the TERAPIX (http://terapix.iap.fr/soft/) center using 
FLIPS and then processed with the TERAPIX modules 
for photometric and astrometric calibrations and mosaic 
image construction. 

In both CRS6S, Oil 'superflat" built from science images 
on "empty" fields obtained during the same period was 
used to subtract the fringe patterns. The mosaic is nor- 
malised on CCD 04, with the higher sky value. Residual 
variations are very small and the flat-fielding appears to 
be good to better than 0.1% in B and V across the full 
mosaic (Kalirai et al. 2001). 

Astrometric and photometric calibration and mosaic 
image construction were performed for each CCD chip us- 
ing IRAF routines for first epoch data. The APM cata- 
log (http://www.ast.cam.ac.uk/apm-cat/) has been taken 
as reference to transform the distorted (x,y) positions 
into corrected (a, (5) coordinates using the IRAF package 
"IMAGE.IMCOORDS" . For second epoch data Terapix 
modules were used. With Terapix, photometric and astro- 
metric calibrations are determined using the entire mosaic 
to give more homogeneous solutions. A global astrometric 
solution for all CCD images projected onto a common sys- 
tem was computed with comparison to the USNO catalog 
for absolute calibration. 

Photometric calibration was performed by comparison 
to a set of several Landolt fields observed at the same 
epoch. For each filter, the zero point determination was 
found to be quite similar for the 12 individual CCDs, with 
a mean error less than 0.025 magnitude, except for the 
South first epoch data in the R filter which will be not 
used in the catalog. 

Because the shifts between exposures are small and 
there are few exposures, gaps between individual CCDs 
are only partially corrected and present a lower signal-to- 
noise ratio. However, this affects only a very small fraction 



of the surveyed area, and should have little or no effect on 
most studies based on the catalog. 

3. The catalog 

3.1. Choice of first/second epoch and North/South 
data 

In the final catalog, data for objects with declinations S < 
-1-28.05° are taken from the South field, while objects with 
S > +28.05° are taken from North data. 

On the whole field, the B and V data are only avail- 
able from first epoch data and the I band data from sec- 
ond epoch data. The R band data are taken from the 
first epoch in the North and from the second epoch in the 
South. 

3.2. Combination of B, V, R and I data 

The BVR catalog for first epoch data was built with ref- 
erence to the R positions, selecting the nearest object in 
the other filters, and checking the magnitude consistency. 
As the astrometric solution was determined CCD by CCD 
using a low order polynomial function, residual distortions 
can be large (1 to 3 arcsec) in the corners and at the edges 
of the individual CCDs, and the proximity criterion is no 
longer valid for combination with data which do not suf- 
fer from the same errors. This is the case for the second 
epoch data for which a better global astrometric solution 
on each of the South and North fields was derived. 

In order to recover the distorted regions, we developed 
a more sophisticated geometrical method, searching for 
the association which minimizes the distances of all ob- 
jects in a given area around the reference objects. This is 
a pattern recognition process which is quasi-independent 
of residual distortions if the area is small compared to the 
distortion scale. In cases when there were fewer than 2 
objects in common in the exploration area, the nearest 
object was selected. The best combinations (in terms of 
errors and number of matched objects), confirmed with 
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a map of the associated positions, was obtained with an 
exploration area of radius 15 arcsec, allowing for position 
errors as large as 3 arcsec. At the end of the process, re- 
sults were checked for possible duplicate associations and 
corrected. 

The final result is a single catalog of more than 60,000 
objects (galaxies and point sources) with B, V, R and I 
data in the area covered by the South and North pointings. 
The object positions were selected from the R data, from 
the first epoch in the North and the second epoch from 
the South. 

3.3. Catalog content 

The following quantities are given in the catalog: 

- identification number for each object 

- {a,5) coordinates in J2000 

- aperture magnitudes in the 4 bands for a 3 arcsec 
aperture 

- SExtractor error estimates on these magnitudes 

- total magnitudes in the 4 bands 

- SExtractor error estimates on these magnitudes 

- SExtractor central surface brightness (in the I and R 
bands) 

- external magnitude errors in the four bands and for 
each measured magnitude 

- shape parameters: minor axis, major axis, orientation 
(in the I and R bands) 

- SExtractor star-galaxy classification (in the I and R 
bands) 

- our star-galaxy classification in the I band (distinct 
from the SExtractor one) 

- SExtractor detection flags (see Bertin & Arnouts 
1996) 

- B, V, and R magnitudes taken from the literature, 
scaled to our system for objects brighter than R=17.5 

How these parameters were computed is described be- 
low. 

4. Catalog analysis 

4.1. Object extraction 

Sources were extracted using the SExtractor package 
(Bertin & Arnouts 1996). The detection threshold was set 
to 2a, a trade off between detecting as many objects as 
possible and limiting the number of false detections. The 
minimum number of contiguous pixels above the detection 
threshold for extraction was set to 9, a conservative value 
which corresponds to about half of the pixels at FWHM 
for a 0.9 arcsec seeing. In this way, although wc limit our 
detection of the large surface brightness population, wc 
increase the confidence in the detections. 

Aperture magnitudes were calculated, in all filters, us- 
ing the SExtractor MAGAuto parameter based on the 
Kron total magnitude. Central surface brightnesses were 
computed in the best quality bands (R and I) only. 



1.5- 




14 15 16 17 18 19 20 21 22 23 24 25 
R maanitude 

Fig. 2. Variation of the astrometric dispersion in a (solid 
lines) and 6 (dashed lines) as a function of R magnitude. 
From bottom to top: Very thin lines (first two lines): com- 
parison between first and second epoch data. Thin lines 
(third and fourth lines): comparison with Bernstein et al. 
(1995). Thick lines (fifth and sixth lines): comparison with 
GP77. 

4.2. Internal astrometry 

The astrometry was computed on the R band images with 
an accuracy of about 0.2-0.3 arcsec over the whole field, 
which is typical of the precision of the reference catalogs, 
except in the small bands at the edges of the CCDs for 
the first epoch data. 

The relative astrometry with the other filters is much 
better; for example, in the B band, the residual dispersion 
is 0.07 arcsec in a and 0.08 arcsec in 6. 

Comparison of the astrometric solutions between the 
two epochs in the common South field and in the R filter 
shows an agreement at the 0.5 arcsec level (see Fig. (S)) for 
all magnitudes. 

4.3. Astrometry comparison to published data 

For sources brighter than R=20, the comparison of our 
catalog positions to those in the USNO and GP77 catalogs 
shows a dispersion of about 0.75 arcsec with USNO and 
0.9 arcsec with GP77 in both a and 5. The astrometry 
dispersion between our catalog and GP77 is typical of the 
astrometry error of the GP77 catalog itself. 

For sources fainter than R=20, the comparison of our 
catalog with that of Bernstein et al. (1995) shows a dis- 
persion of about 0.7 arcsec. 

Results are presented in Fig. |21 

4.4. Photometric corrections 

Zero points were computed using observations of stan- 
dard stars (Landolt 1992) in the fields SAlOl, SA104 and 
SAllO. These standards were observed at nearly the same 
air mass as the science observations in the case of first 
epoch data and scaled to zero air mass for the second 
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epoch. The correction in magnitude due to small airmass 
variations would have been, at maximum, of 0.033 magni- 
tude for R, 0.020 magnitude for V, and 0.015 magnitude 
for B, much less than the estimated error in the magni- 
tudes, and hence were neglected. We also computed cen- 
tral surface brightnesses in the R and I bands. The cali- 
bration was done on the entire mosaic after scaling each 
chip, assuming identical color equations for all CCDs. This 
proves to be true at better than 4% for all the CFH12K 
observations (McCracken et al. 2003). The B, V and R 
filters appear to have negligible color terms with respect 
to the Johnson Kron Cousins system and we chose to stay 
in the CFH12K system, applying no correction in the cat- 
alog for the color equations, even in I where the correction 
could be of the order of 1=0.1. Instrumental magnitudes 
are converted to the Vega system. 

We correct for Galactic extinction using the Schlegel et 
al. (1998) maps. This contribution remains small: ~0.05 
magnitude for the B band in the worst case. 

4.5. Magnitude error budget 

Magnitude errors are due to zero point uncertainties, see- 
ing inhoniogeneities and internal errors. These are dis- 
cussed below. 

4.5.1. Computed error budget 

Zero point uncertainties are smaller than 0.03 magnitude. 
In order to check the effects of seeing conditions on the 
photometry, we compared the aperture magnitudes on R 
images (0.8 arcsec seeing) to the aperture magnitudes on 
the same images convolved with a Gaussian to mimic 
1.00 arcsec seeing. We confirm the results by Saglia et 
al (1993): when the seeing gets worse, the objects ap- 
pear fainter since part of the signal is scattered into the 
background. The worst error remains smaller than 0.25 
magnitudes at R^24. We note, however, that such a pro- 
cess slightly underestimates the error because of the back- 
ground smoothing. 

Internal SExtractor errors are of the order 0.1 magni- 
tude at R~24, a large part probably due to close neighbors 
as shown in the section which describes the catalog prop- 
erties. 

In total, we therefore expect errors for the faintest cat- 
alogued objects (R'--^24) to be of the order of ^0.27 mag- 
nitude. 

4.5.2. Externally estimated error budget 

To get an external estimate of the magnitude error budget, 
we used the overlap between the North and South fields in 
the B, V and I bands and the two South pointings in the R 
filter. Results based on several thousand objects common 
to both fields in the overlap region for the B and V filters 
are shown in Figs. Eland 0as a function of magnitude. 




Fig. 3. Upper graph: solid line: statistical la magnitude 
uncertainties between the North and South fields, dashed 
line: quadratic sum of the internal error from SExtractor 
and of the error due to the different seeings between dif- 
ferent observations. Lower graph: B North versus B South 
magnitudes. 




Fig. 4. Upper graph: solid line: statistical Icr uncertainties 
between North and South fields, dashed line: quadratic 
sum of internal error from SExtractor and error due to 
the different seeing between different observations. Lower 
graph: V North versus V South magnitudes. 

Results of the two R observations and in the common I 
area are shown in Fig. [3 

The increasing uncertainty at the bright end in the 
R and I bands is due to the second epoch data, because 
of a background over-correction by Terapix near bright 
objects, and saturation effects on deep images. 

^From this analysis we can approximate the magnitude 
dispersion by regression laws for B, V, R and I filters with 
parameters given in Table [S] 

The photometric accuracy is better than 0.3 magni- 
tude in all bands, in good agreement with the computed 
error budget. Also, the R error estimate is probably an 
overestimate for the faint part of second epoch data, which 
are clearly deeper and of better quality. 
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Fig. 5. Upper figure: Statistical la uncertainty between 
first epoch and second epoch data on South R field. Lower 
graph: Statistical Icr uncertainty between I band second 
epoch data in the common area. 

Table 5. Fit of the magnitude uncertainty as a func- 
tion of magnitude deduced from the common areas. This 
was used to compute uncertainties for individual galaxies: 
Uncertainty = slope x magnitude + constant term. 



Magnitude band 


validity domain 


slope 


constant term 


B 


[20,25.5] 


0.0340 


-0.5933 


V 


[20,24.5] 


0.0445 


-0.8352 


R 


[20,24] 


0.0093 


-0.0160 


I 


[20,23.5] 


0.0657 


-1.1772 



4.6. Magnitude biases 

Systematic bias in the magnitudes of catalogued objects 
could occur as a result of the presence of close neighbors 
and saturation effects. The effects of each of these is as- 
sessed below. 



4.6.1. Neighbor contamination 

To check this effect on the magnitude determination, we 
select the objects flagged by SExtractor as having a chance 
of being biased by at least 0.1 magnitude by a close neigh- 
bor. The percentage of such objects, measured in the R 



band, is shown in Fig. Eland appears to be at a relatively 
constant level close to 15-20%. This tells us that a rela- 
tively small fraction of our objects is significantly affected 
by close neighbors. 

4.6.2. Saturation effects 

In order to quantify the percentage of objects which could 
be affected by saturation effects, we select objects flagged 
by SExtractor as having at least one pixel saturated in 
their profile. The result is given in Fig. as a function of 
magnitude and shows that up to R=:17.5 objects have a 
probability to have at least one saturated pixel higher than 
50% . The magnitude error is not quantified and depends 
on the object brightness and morphology and so we chose 
to check the magnitude of all brighter objects in shallower 
published data. 

We have added published data only if they cover our 
whole field of view (GP77, Terlevich 2001 and Komiyama 
2002) to our catalog of objects brighter than R=17.5 ap- 
plying the mean magnitude shifts computed in Section 4.7, 
without color corrections. 

In the V band data, objects from the Terlevich et al. 

(2001) catalog fainter than V~14 should not be saturated 
and can therefore cover the whole range of our saturated 
objects. 

In the B and R bands, data acquired by Komiyama et 
al. (2002) are supposed to not be affected by saturation 
for magnitudes fainter than R=14. However, we expect 
that very concentrated objects could still be saturated for 
magnitudes somewhat fainter than R=14, and so we chose 
a conservative approach by using the Komiyama et al. 

(2002) catalog only between R=15.75 and 17.5 and the 
GP77 catalog for brighter objects. 

In summary, all objects in the catalog fainter than 
R=17.5 come from our CFH12k data. For objects be- 
tween R=15.75 and 17.5, their R and B magnitudes are 
taken from Komiyama et al. (2002) or GP77 if the objects 
are classified as saturated by us. Objects brighter than 
R=15.75 are taken from GP77 if classified as saturated 
in our data. Finally, V magnitudes of objects brighter 
than R=17.5 and classified as saturated are taken from 
Terlevich et al. (2001). 

In total we identified ~540 potentially saturated ob- 
jects in our catalog with R<17.5, of which 272 are iden- 
tified as galaxies in the literature. The unidentified ob- 
jects could be either stars (which are not included in the 
Komiyama et al. and GP77 catalogs) or galaxies. However, 
at least half of the R<17.5 objects are stars (see Fig. I15|) . 
which are more easily saturated than galaxies, and hence 
the largest part of the unidentified objects is likely to be 
stars. 

We looked, however, at the 270 unidentified objects in 
our images and most of them were obviously stars (slightly 
less than 70%). A smaller subsample was made of com- 
pact objects, possibly stars or spheroidal galaxies (slightly 
less than 30%). Finally, a few galaxies were identified (7 
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14 15 16 17 IB 19 20 21 22 23 24 
R 



Fig. 6. Solid line: percentage of objects with at least one 
saturated pixel as a function of R magnitude. Dashed line: 
same for biased magnitudes by 0.1. 
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Fig. 7. Solid line: statistical la error between R GP77 
and R CFHT magnitudes, dashed line: B GP77 versus B 
CFHT magnitudes. 
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Fig. 8. Statistical Icr uncertainties between V magnitudes 
from Terlevich et al. and from our CFHT observations. 
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Fig. 9. Solid line: statistical Icr uncertainties between R- 
band magnitudes from Komiyama et al. and our CFHT 
data, dashed line: statistical Icr uncertainties between 
Komiyama et al. and CFHT B magnitudes. 



galaxies). We therefore decided to classify as stars all the 
unidentified objects that were not obviously galaxies. This 
number of ~263 stars is, moreover, in good agreement with 
the Besan^on model star counts (Gazelle et al. 1995). 

4.7. Magnitude comparison with other catalogs 

We first compared our catalog to the GP77 photographic 
catalog (up to R~19.5) which covers a very large area that 
encompasses our whole field. The result is given in Fig. [3 
The dispersion is lower than 0.4 magnitude at R=20 and 
close to 0.1 in R and 0.15 magnitude in B at R=15.5, with 
the following mean offsets: 

B CFHT = B GP77 + 0.20 

R CFHT = R GP77 + 0.52 

We also compared our results with those of Terlevich 
et al. (2001) in the V band and Komiyama ct al. (2002) 
in the B and R bands. The results are shown in Figs. |H1 
and El B, V and R data are in good agreement, with a 



dispersion at the bright end close to 0.2 magnitude for V, 
0.25 magnitude for R and 0.3 magnitude for B. The mean 
offsets are; 

B CFHT = B Komiyama - 0.09 

V CFHT = V Terlevich + 0.12 

R CFHT = R Komiyama + 0.03 

Finally, we compared our results with the CCD ob- 
servations of Bernstein et al. (1995). The result is shown 
in Fig. ^1 At R=20, the dispersion is within 0.15 mag 
for the R band and 0.2 mag for the B band. At the faint 
end of the catalog (R^25) the dispersion is 0.35 in R and 
0.45 in B. At our catalog completeness level (R~24) the 
dispersion is 0.25 in R and 0.3 in B. We note that there 
is a moderate to strong difference in seeing between our 
data and those of Bernstein et al., since the worst seeing 
for our data was 1.07 arcsec in B and 0.87 arcsec in R, 
while that of Bernstein ct al. was 1.31 arcsec. This proba- 
bly explains why we see a rather smooth variation of the 
dispersion between the Bernstein et al. and our B band 
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Fig. 10. Solid line: statistical la uncertainties between 
R Bernstein et al. and R CFHT magnitudes, dashed line: 
statistical la uncertainties between B Bernstein et al. and 
B CFHT magnitudes. 

magnitudes, while the comparison in the R band shows 
visible oscillations. 

The consistency between our catalog and others allows 
us to merge some of the data to recover the saturated 
bright end of our catalog. As we will use external data 
only for bright galaxies over a rather limited magnitude 
range, and the mean offsets are small, we do not include 
color terms from one photometric system to an other. 

5. Catalog completeness 

The catalog completeness has been measured in two ways: 
by simulations and by comparison with a deeper catalog. 

The simulation method adds artificial objects of differ- 
ent shapes and magnitudes to the CCD images and then 
attempts to recover them by running SExtractor again 
with the same parameters used for object detection and 
classification on the original images, fn this way, the com- 
pleteness is measured on the original images and at differ- 
ent locations in the cluster. We investigated the catalog 
completeness for point-like and faint low surface bright- 
ness objects separately. This is because part of the faint 
galaxy population in clusters consists of faint low surface 
brightness galaxies (R central surface brightness fainter 
than 24), which can be nucleated (e.g. Ulmer et al. f996). 
These objects are crucial for understanding the cluster 
physics and, determining how deep we can observe this 
population by comparison with point like objects, is of 
major interest. 

5.1. Simulations for point-like objects 

The results at the 90% mean completeness level are sum- 
marized in column two of Table El for all filters and in the 
two fields. 

An example is given in Fig. If f I which shows how the 
completeness levels vary from CCD to CCD because of the 



QE variations between individual CCDs and because of 
fluctuations in the diffuse background light due to bright 
stars and galaxies in the field. 

5.2. Simulations for low surface brightness objects 

We estimated the completeness of our catalog for low 
surface brightness galaxies using simulated point-like ob- 
jects with a FWHM of slightly more than 2 arcsec from 
a Gaussian profile. This is the typical maximal size of a 
low surface brightness galaxy (e.g. Ulmer et al. f996) in 
Coma. This method does not always take into account the 
true profile of low surface brightness objects in Coma but 
is a good compromise between simulation simplicity and 
result accuracy. Results are summarized in column 3 of 
TableEland an example is given in Fig.^lfor the B filter. 

5.3. Comparison with the Bernstein catalog 

Assuming that the Bernstein catalog is fOO% complete 
at our limiting magnitude, we can use it to obtain an 
independent estimate of the completeness of our catalog. 
Based on this, we find 90% completeness at R=23.5, which 
differs from that predicted by our simulations, 90% at 24.5 
for point-like objects and 50% at 22.5 for faint low-surface- 
brightness objects. 

There are probably several reasons for the discrepancy. 
First, faint galaxies are not exactly stellar-like objects and 
their detection rate is, therefore, lower than for stars. This 
probably explains part of the difference. 

Second, the Bernstein (f 995) field is located very close 
to the two dominant galaxies of the Coma cluster where 
there is a high level of diffuse light, which can affect the 
detection of faint objects against this background. The 
Bernstein image was specifically treated to remove this 
light and improve the detection rate of the faintest objects. 
As we do not try to correct our data from diffuse light, this 
may contribute to the difference in completeness limits. 

6. Star-galaxy separation 

Objects are classified as star-like or galaxies following two 
criteria both based on the f band data. We used this band 
because it was reduced in the most homogeneous way on 
the whole field and was obtained under the best seeing 
conditions. The price to pay is that not all objects detected 
in R have an I counterpart (R band data are deeper than 
I band data). However, as we will limit the star-galaxy 
separation to I=2f , this is not a serious concern as most 
of the l<2f objects are detected in the 4 bands. 

As a first criterion for star-galaxy separation, we use 
the SExtractor class parameter, which varies from for 
galaxies to f for stars. We classify as stars objects with 
values greater than 0.98 (e.g. McCracken et al. 2003). 

The second criterion comes from the relation between 
the total magnitude and the central surface brightness for 
different types of objects. The stellar locus clearly appears 
up to I~21 on Figs. Eland im showing unambiguously 
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the good quality of the star-galaxy separation up to this 
magnitude. The lines separation between stars and galax- 
ies were put in by hand in order to optimally separate the 
stellar and galaxy loci. 

For objects between 1=17 and 18.5, we did not use the 
total magnitude / central surface brightness criterion be- 
cause the separation line starts to classify stars as galaxies 
around 1=18.25 and classify as stars all objects based on 
the SExtractor class criterion. For objects between 1=18.5 
and 21 we used the second criterion and classified as stars 
all objects below the heavy solid line in Fig. E| 

Because the North and South data were obtained un- 
der different seeing conditions, the criteria were slightly 
different in the two fields and so we performed the star- 
galaxy separation individually for each of the fields. 

Beyond 1=21, the star-galaxy separation is very un- 
reliable, because of possible confusion of stars with small 
galaxies (with seeing dominated profiles) and we chose to 
classify all objects in this range as galaxies. 

6.1. Star counts: comparison to the Besanfon model 

Our star counts shown in Fig. 1151 agree well with the pre- 
dictions from the Besangon model for our galaxy (Gazelle 
et al. 1995 or http://bison.obs-besancon.fr/modele/) in 
the I band and confirm the quality of the star-galaxy sep- 
aration between 1=17 and 1=21. At fainter magnitudes, 
the predictions show that the contribution of stars is only 
10% at 1=21 and drops quickly beyond that. 

6.2. Moving objects and spectroscopically confirmed 
stars 

In order to further check the position of the stellar locus in 
the surface brightness/total magnitude plots, it is useful 
to plot known stars in the diagram. One way to do this is 
to use proper motion to identify stars in the Galactic halo. 
However, a check of the Vizier database at CDS revealed 
that there are no published moving stars in the I=[19;21] 
range in this field. 

We therefore compared our images observed in 2000 
to the DSS photographic plate observed in 1955 and to 
the Bernstein et al. (1995) field observed in 1991, in or- 
der to identify moving stars in our fields. Using a sim- 
ple cross-correlation we select all CFHT objects with a 
counterpart between 3.4 arcsec (i.e. ~3 times the rela- 



tive astrometric precision between catalogs) and 15 arc- 
sec in the DSS and between 1.4 arcsec (i.e. 2 times the 
relative astrometric precision between the catalogs) and 
4 arcsec in the Bernstein et al catalog. These upper val- 
ues limit our search to proper motions smaller than 400 
mas/year, a reasonable upper limit (see e.g. the Tycho-2 
catalog. Hog et al. 2000). After examining each of these 
objects and rejecting spurious detections, we were left 
with two bona fide moving objects in the South field and 
none in the North field. The first of these is detected 
from the DSS plate (a2ooo = 194.8629°, (520oo=27.6834°, 
1=19.62). The second one is detected from the Bernstein 
et al. (1995) image (a2ooo = 194.8673°, 52ooo=27.8909°, 
1=19.67). Additionally, observations reported in Adami et 
al. (2000) led to the spectroscopic discovery of a third star 
(a2ooo=194.8149°, (52ooo=27.9257°, 1=19.05) in our fields. 

These three objects are well located in the star locus, 
giving us confidence in the reliability of our star-galaxy 
separation at least down to I~20. 

7. Colors 

In order to check the reliability of the photometric spec- 
tral energy distributions (SEDs) obtained in our fields, we 
have compared the observed colors with theoretical ex- 
pectations derived for two different and well-defined sam- 
ples of objects: stars and galaxies with known redshifts. 
Although photometric redshifts and star-galaxy identifi- 
cations are beyond the scope of the present paper, this 
section provides a reference on the quality of our SED 
data from ~ 4000 A to 1 /xm. 

7.1. Stars 

Synthetic colors for stars have been derived for a variety 
of spectral types and luminosity classes, using the empir- 
ical stellar library of Pickles (1998). We have adopted a 
detailed modeling for filter transmissions, taking into ac- 
count the total efficiency of the system as a function of 
wavelength. Stars have been selected in our catalogs ac- 
cording to the criteria given in Sect. 7 between 7=17 and 
21. There is a good agreement between synthethic and ob- 
served colors of stars, as shown in Fig. ^] On the color- 
color diagram B-V vs. R-I, both the theoretical locus of 
the observed main sequence and the expected dispersion 
towards B-V > 1.2 are well reproduced by observed stars. 
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Fig. 11. Completeness in percentages in B, V, R and I 
magnitudes for point-like objects in the North field. The 
dotted lines show the completeness of individual CCDs, 
the solid line is the mean completeness for the North field 
and the solid horizontal line is the 90% completeness level. 
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Fig. 12. Same as Fig. ^2 for the faint low surface bright- 
ness galaxies in the North field and in the B band. 

Fig. 13. Central I surface brightness versus I total magni- 
tude. Crosses are objects identified as stars by SExtractor 
with class parameter 0.98 or greater. The heavy solid line 
shows the star-galaxy separation we chose for the range 
I=[18.5;21]. Upper graph: North field. Lower graph: South 
field. 

Fig. 14. Same as Fig. [T31 but in I=[18.8;22]. The large 
squares in lower graph are the 3 faintest bona fide stars 
(determined by external means) in the field. 




17 IB 19 20 21 22 

I 



Fig. 15. Star counts in our Coma catalog (long-dashed 
line), star counts from the Besangon model (thin line with 
error bars) and galaxy counts (solid line). Error bars are 
at the 3tT level. 

7.2. Galaxies 

The set of empirical templates compiled by Coleman, Wu 
and Weedman (1980) was used to derive representative 
synthetic colors for galaxies in our photometric system as a 
function of redshift. The spectro-morphological types con- 
sidered are the following: E, Sbc, Scd and Im, thus provid- 
ing a simple representation of the galaxy population in the 



Adami et al.: Deep and wide field imaging of the Coma cluster: the data 



11 




1 2 

B-'J 



Fig. 16. R-I versus B-V for I=[17;21] stars in our catalog 
(dots), compared to synthethic values derived from the 
library of Pickles (1998) (large open circles). 



local universe. Colors were obtained using the detailed fil- 
ter transmissions mentioned above for the CFH12K cam- 
era, with a simple ^-correction. 

A sample of 873 cluster members compiled from the lit- 
erature (with measured velocities between 4000 and 10000 
km/s) is available, as described by Adami et al. (2005b). 
Although this spectroscopic sample is rather small, ob- 
served colors are found to be in good agreement with 
model predictions. A representative example is given in 
Fig- El for B-I as a function of redshift. A large major- 
ity of cluster galaxies display colors fully compatible with 
early-type models, and the range of colors spanned by the 
whole spectroscopic sample is in agreement with simple 
model expectations. This is particularly true for galaxies 
morphologically identified as early types from visual in- 
spection (Biviano et al. 1996). It is worth noting that dust 
reddening has not been considered here, and some objects 
displaying extremely red colors would need this correction 
to fit into the scheme. In particular, a moderate intrisic 
extinction Ay > 0.5 magnitudes (E(B-V) > 0.12, with a 
Calzetti et al. 2000 reddening law) provides B-V > 2.5 
when applied to an Sbc template at the Coma redshift, 
using a simple dust-screen model. Thus, dust extinction 
could naturally explain the very red colors observed for 
a few galaxies (less than 10% of the total spectroscopic 
sample) in Fig. El compared to models. Also bright ob- 
jects with pixels close to the saturation limit could exhibit 
atypically red colors. 

Also bright objects with pixels close to the saturation 
limit could exhibit atypically red colors. 

In summary, the bulk of the spectroscopic sample 
nicely fits into the expected color-redshift diagrams, and 
this indicates a good quality of the photometric SEDs for 
future studies based on SED-fitting techniques. 




0.2 0.4 0.6 o.a 




1 Li , I d 

0. 2 

z 

Fig. 17. B-I versus redshift for a sample of galaxies with 
spectroscopic redshifts in this field (from Adami et al. 
2005b), compared to simple model expectations derived 
from Coleman et al. (1980) templates. From top to bot- 
tom: elliptical galaxy (solid line), Sbc (dotted Hne), Scd 
(long dashed line) and Im (dot-dashed Hne). Top panel: 
all galaxies. Bottom panel: galaxies visually classified as 
early type galaxies according to their morphology. 

7.3. Red sequence in the Color Magnitude Relation 

We also used the spectroscopic sample described in the 
previous section to investigate the color magnitude re- 



Fig. 18. Color magnitude relation in the Coma cluster. 
Small dots: all galaxies on the Coma line of sight. Grey 
(B&W format) / green (color format) filled circles: Coma 
members. Black (B&W format) / red (color format) filled 
circles: early type Coma members. Black (B&W format) / 
red (color format) solid line: our fitted red sequence. Upper 
figure: grey (B&W format) / green (color format) solid 
and dashed lines are the 1 and 3cr error envelopes. Lower 
figure: on an expanded scale, the solid grey (B&W format) 
/ green (color format) line is the la error envelope and 
the two black (B&W format) / red (color format) dashed 
lines are the red sequences from Seeker et al. (1997) and 
Lopez-Cruz et al. (2004) 
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lation (B-R versus R, hereafter the CMR) in the Coma 
cluster and to compare our results with published red se- 
quences in Coma. The Coma cluster is known to exhibit 
a well defined red sequence (e.g. GP77 or Mazure et al. 
1988). The shape of this relation is probably driven by 
metallicity effects in different mass systems (e.g. Kodama 
& Arimoto 1997). Massive galaxies can be redder than 
lower mass systems because they are able to retain more 
metals and then to form more evolved stars. Early type 
galaxies in this quite old cluster (e.g. Adami et al. 2005b) 
therefore form a well defined sequence with a negative 
slope (faint and low mass early type galaxies retain fewer 
metals than bright and massive ones and are therefore 
bluer) . 

Fig. 1 181 shows this relation. We fit a red sequence using 
only galaxies classified as early type galaxies and likely 
Coma members based on their velocities (between 3500 
and 10000 km/s). The relation is: 

B-R = -(0.045±0.028) x R + (2.27±0.48) 
The errors are quite large, but this is simply due to 
the small size of the spectroscopic sample (even if one of 
the largest available for a cluster of galaxies: Adami et al. 
2005b) and to the uncertainties in the visual morpholog- 
ical classification (potentially including late type galaxies 
in the sample). We also overplot on Fig. ^| the error en- 
velopes at the 1 and 3a levels around this mean relation. 
Error envelopes were simply computed using the modeled 
magnitude uncertainties (see Section 4.5) as a function of 
magnitude. 

Several points are noteworthy. First, the red sequence 
defined in this way is in good agreement with the galaxy 
distribution in the R/B-R space. The high galaxy den- 
sity regions in this diagram are well correlated with the 
red sequence. We also visually inspected the two early 
type galaxies inside the Coma cluster with extreme B-R 
values. These two objects (around B-R=2.5, R=16.3 and 
B-R=0.8, R=18.3) have several velocity measurements in 
the literature (5 and 4 respectively), all consistent with 
them being Coma members (7451 km/s and 8043 km/s). 
These objects must certainly have had peculiar histories 
that made their B-R colors atypical. They probably ex- 
perienced dynamical encounters at least for the reddest 
one which has a peculiar shape in our images. However, 
discussing the precise evolution of these two objects is be- 
yond the scope of this paper. 

Second, the galaxies inside the Coma cluster but with- 
out morphological information available in the literature 
(filled grey circles in Fig. I18|) are also very close to this 
red sequence, defining a relatively narrow relation down 
to R~19 (consistently with e.g. Adami et al. 2000). What 
occurs at fainter magnitudes remains uncertain as we do 
not have deep enough spectroscopy to reach any conclu- 
sion (see also Adami et al. 2000). 

Finally, our red sequence agrees very well with results 
published previously in the literature, as shown in Table 
[7| Some marginal differences exist for the constant terms 
but are due to slight differences in the B and R filters used 
by different studies. 



Table 7. Literature red sequences in the B-R/R sequence. 



Authors 


Slope 


Constant term 


present paper 


-0.045 


2.27 


Gladders et al. (1998) 


-0.045 




Lopez-Cruz et al. (2004) 


-0.046 


2.22 


Seeker et al. (1997) 


-0.056 


2.41 



Investigations of the faint part of the red sequence are 
beyond the scope of this paper and will be addressed in a 
future work. 

8. Summary and prospects 

We have presented and discussed the properties of a cata- 
log of more than 60,000 objects in the Coma cluster based 
on deep CFHT observations in the B, V, R and I bands 
and spanning an unprecedented 10 magnitude range in one 
of the largest ~0.72x '^0.82 deg^ areas presently available 
for Coma at this depth. The catalog is complete at the 90% 
level down to R'-~^24 for stellar-like objects and to R~20.75 
for faint low surface brightness (see Sect. 5.2) galaxy-like 
objects. 

Astrometry is accurate to 0.5 arcsec, except for first 
epoch data on the edges and corners of the CCDs. 
Magnitude errors are smaller than ~0.3 in all bands 
down to R=25. Saturation effects on objects brighter 
than 1=17.5 have been removed by using shallower pub- 
lished data. More generally, galaxy magnitudes are in good 
agreement with published data and the colors of galaxies 
and stars are in good agreement with synthetic models. 
The bright part of the CMR agrees well with that derived 
for galaxies with known redshifts in the Coma cluster, and 
with previous CMRs published in the literature. 

The star-galaxy separation is robust for all objects 
brighter than 1=21, and the star counts fit the Besangon 
model very well. 

These data have already been used to search for dif- 
fuse emission in the Coma cluster (Adami et al. 2005a) 
and to look for and analyze properties of faint low surface 
brightness galaxies (Adami et al. 2006). We plan next to 
use these data to investigate the properties of the different 
cluster galaxy classes and to derive the luminosity func- 
tions of the Coma cluster galaxies in various bands and 
in various regions of the cluster. This should deepen our 
knowledge of environmental effects on galaxy luminosity 
functions, already shown to be strong in Coma (e.g. Lobo 
et al. 1997). 

The catalog described above will be made public by the 
CENCOS center at http://cencosw.oamp.fr/. Individual 
catalogs for the second epoch R and I band data will also 
be available upon request. 
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